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The primary purpose of this thesis was to provide a 
method of checking the geometry and element connectivity 
input data for two finite element programs, ADINA and 
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r.  LNTRODUCTION 


A. GENERAL 

The continued developments and advancements of the 
finite element method this last decade have provided 
greater machine capabilities than ever before. Here at 
the Naval Postgraduate School (NPS), two of the more 
favorable finite element programs are SAP IV [reference 1] 
for linear analysis and ADINA [reference 2] for nonlinear 
analysis. However, with the large amounts of numerical 
input/output data and automatic mesh generation, it is 
impractical to check and reduce this data without a graphi- 
cal representation. Data checking is divided into two 
areas: preprocessing and postprocessing. 

i. Preprocessing 

Preprocessing is the checking of the input data 

deck. Errors in a finite element program occur basically 
in two areas. First, how close is the mathematical model 
(boundary conditions, loading conditions, material proper- 
ties, etc.) to the real problem? Second, are numerical 
errors present, or did misinterpretation of instructions 
occur in data deck preparation? Of the second type, most 
common errors are found in the geometry and element connec- 


tivity data. Preprocessing includes the forming of a 


which geometry and element errors are easily detected. When 
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node and element numbering options are available, the graph 
aids in the physical interpretation of the output data. 
Preprocessing is not a foolproof method of eliminating 
errors, but it does provide a tremendous advantage to the 
user. Preprocessors may be incorporated into the data check 
mode of the finite element program. However, to modify a 
large and complex program is dangerous. This may not be 
the best approach. A safer method is to develop a pre- 
processor which will read the finite element program deck 
separately with minimum modifications to that deck. 
2. Postprocessing 

Though not as important as preprocessing, post- 
processing is extremely helpful in output analysis. Probably 
the most common and useful type of postprocessing is the 
contour plot. Appendix C of reference 4 lists a program 
developed to produce contour plots of stress data from finite 
element models. Contour plots can easily be adapted to a 
2D system, but 3D requires plotting the contours on 2D 
surfaces, a bit more complicated. Two methods of post- 
processing of displacements are the plotting of a deformed 
model or placing scaled vectors at the nodes. In the cases 
where the displacements are small, multiplication by a 
magnification factor produces an exaggerated representation. 
Like preprocessing, postprocessing can be incorporated into 
the finite element program directly, or done separately. 
When done separately, the finite element program must still 


be modified slightly to obtain a punched deck of the stresses 
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and displacements in the desired format. This thesis is 


primarily concerned with preprocessing. 


B. HISTORY OF DEVELOPMENT OF PSAP1 

A package [reference 4] containing digital computer 
programs for generating oblique orthographic projections 
and contour plots was produced by the National Aeronautics 
and Space Administration's (NASA) Langley Research Center 
(LRC) and distributed by the National Technical Information 
Service in January, 1975. The programs are completely 
general. Both programs contain options for selecting 
various plotting equipment including CALCOMP, VARIAN, and 
cathode ray tube (CRT) displays. With minor modifications, 
they can be adapted to most any system. 

1. SUBROUTINE PSAP Implementation 

Losh [reference 6], for his master's thesis in 

aeronautical engineering, implemented the preprocessor and 
postprocessor program, PSAP, at NPS in December, 1976. 
Modifying the LRC package [reference 4], Losh adapted 
SUBROUTINE PSAP to the NPS IBM 360/67 system using the 
CALCOMP model 765 plotter. PSAP serves as a preprocessor 
for SAP IV models, and serves as a postprocessor for dis- 
placements of those models. Unfortunately, PSAP is severely 
limited in the type of elements it can plot. 

2. Motivation for SUBROUTINE PSAPI1 

With the introduction of the ADINA [reference 2] 


program at NPS in January, 1977, and with expectation of 
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doing future analysis on ceramic turbine blades, it was 
desired to expand PSAP to include all ADINA elements and 
the 8-20 node brick elements in SAP IV. Like PSAP, PSAP1 
contains preprocessing and displacement postprocessing 
capabilities. PSAP1l is presented in this thesis and has 
the following improvements over PSAP: 

a. Preprocessing for all ADINA elements. 

b. SAP IV 8 and 8-20 node elements. 

c. Expansion of SUBROUTINE ERROR. 

d. Interpolation of curves using shape functions 
[reference 3] through three points on the edges of the 
8-20 node brick elements and the 4-8 node plane elements. 

e. Improvements in defining the plot origin. 

ESNEAGEZELONOOLE anioptien’ (ISCALE = 0) to plot 
sections of a model without losing perspective. 

g. Several other minor modifications considered 


improvements. 


C. . PRESENT CAPABILITY 

PSAP1 has the capability to plot all ADINA elements 
and all SAP IV elements except the pipe element. It will 
interpolate curves on the edges of 4-8 node plane elements 
and 8-20 node brick elements. Many options are given in 
Appendix A. Some of the more frequently used options are 
listed below: 

1. Numbering of grid points (NOTAT = 1). 

2. Numbering of the elements (NOTAT = 2). 


33 Exploded plot (KDISP = 2). 
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4.  Postprocessing of displacements (NUDISP or NVDISP 
or NWDISP = 1) in two forms: plot of deformed structure 
(KDISP = 1) or displacements represented by vectors at 
the nodes (KDISP = 3; see reference 6). 

5. Symmetric representation about the XY (KSYMXY = 1), 
XZ (KSYMXZ = 1) or YZ(KSYMYZ = 1) planes. 

6. Option to plot sections of the model (partial plot) 
to obtain a better view. Partial plots may be plotted to 
the scale of the complete model to avoid losing perspective 
(ISCALE = 0) or blown up to obtain a better view (ISCALE = 1). 
Multiple plots may be obtained using the same geometry and 
same displacement data (KODE = 1), same geometry and new 
displacement data (KODE = 2), or new geometry and new 
displacement data (KODE = 3; see figure l). 

In general, multiple plots (sections, partial plots, 
additional problems) present no problem. Plotting package 
user courtesy dictates that no more than 5 plots be placed 
on the CALCOMP plotter at any one time. Also, if the plots 
contain many elements (especially 8-20 node elements), it 
is possible to run out of space in the plotting data sets. 


When this happens, you will receive 
ERROR IHC240I STAE, ABEND CODE IS: SYSTEM OB37 SYSPLOT. 


The best thing to do is split the run into two jobs. If 
the job must be run on one job (i.e., a large number of 
elements in the model or an assembly drawing where the 


scale of multiple plots is the same), then SYSPLOT space 
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may be increased [references 7 and 8] by adding the card 
//GO.SYSPLOT DD SPACE=(CYL, (needed space) ) ,SYSOUT=C 
qust.prior.to. card 

//GO.FT1LOFOO1 DD UNIT = SYSDA 


in Appendix A. It would be wise to seek advice from a 
consultant in Ingersoll 146 if additional plotting space 


is required. 


D. EASE OF MODIFICATION 

Both PSAP and PSAP1 are written to maintain as much 
generality as possible for ease of expansion and modifica- 
tion. Several FORTRAN statements, variables and subroutines 
are not used. They were left purposely unchanged. Although 
PSAP1 specifically reads ADINA and SAP IV data, it can 
easily be expanded to include any geometry and element data 
format. Simply study the read-in and storage methods (see 
Section II), and construct appropriate subroutines to read 


amy particular format. 
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II. PROGRAM ORGANIZATION AND DESCRIPTION OF OPERATION 


A. PSAP1 FLOW CHART 

Figure 1 is a condensed flow chart of PSAP1. Probably 
the most important information given on this chart is the 
sequence in which the data cards, NAMELIST OPTION and 
NAMELIST PICT are read. Remember, when generating a 
sequence of plots, once a parameter has been defined, it 
retains that value until it is reassigned. Note that when 
KODE = 1 or 2, the original values of NAMELIST OPTION and 
NAMELIST PICT are retained until they are changed. However, 
when KODE = 3 a new title card, NAMELIST OPTION and a set 
of problem data are read. All variables in NAMELIST OPTION 
and NAMELIST PICT are assigned their default values. A 
new problem begins in this case. It is important that the 
last NAMELIST PICT to be read must contain the value of 


KODE = 0. 


B. NAMELISTS AND EULER ANGLES 
1. NAMELIST OPTION 

Description and default values of NAMELIST OPTION 
are given in Appendix A. Basically NAMELIST OPTION varia- 
bles pertain to the given problem: the number of nodes, 
geometry format, displacement format, space between plots, 
and paper size. Default values for NAMELIST OPTION are 
set, and NAMELIST OPTION is read at the beginning of the 


problem. Assigned values will remain until exit from the 
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Figure 1. Flow chart 
of program PSAP1 


Open plotting data sets. 
Define the origin. 


Read and plot title card. 

Assign defaults to s&OPTION 
and &PICT. 

Read &OPTION. 

Read and print geometry 

and element data. 


KUDISP=KVDISP= 


Read and print displace- 
ment data. 


Rotate original coordinates. 
Plot the elements. 


Plot tvele at top of graph. 


ODE=0 
Position the pen. 

Call PSTOP and close data sets. 
Stop 
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program (KODE = 0) or a new set of problem data is read 
(KODE = 3). 
2.4 NAMELIST PICT 

Like NAMELIST OPTION, NAMELIST PICT variable des- 
criptions and default values are given in Pena A. 
Basically NAMELIST PICT variables pertain to a given plot. 
One NAMELIST OPTION may apply to several successive plots, 
but each NAMELIST PICT defines a unique plot. That plot 
may include the whole model, part of the model and any 
options defined in NAMELIST PICT. NAMELIST PICT also speci- 
fies the viewing plane through the Euler angles (figure 18, 
Appendix A). 

a. Oblique Orthographic Projections (Euler Angles) 

An example of an oblique orthographic projection 

of a finite element model is given in figure 18 in Appendix 
A. The model can be viewed in any selected orientation. 
Euler angle transformations are used to specify orientation 
of the model to be projected. As described in reference 4, 
this transformation resolves the coordinate system of the 
model to a principal viewing plane (i.e., 2 ee Sie) 
em which the display is to be plotted. Prior to rotation, 
the model coordinate system (X,Y,Z) is coincident with the 
coordinate system containing the viewing planes (Xo7% 5725) - 
The viewing planes are fixed, and the model is rotated about 
its model coordinate system. The rotations (,6,6) of the 
body about the model axes (X,Y,Z) are shown in figure 18, 


Appendix A. The NAMELIST PICT variables KHORZ (horizontal 
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axis), KVERT (vertical axis), PSI (ww), THETA (6), and 
PHI (9) specify the viewing plane and Euler angles. The 
order of the Euler angle rotations must be PSI, THETA 


and then PHI. Mathematical transformations are: 


X X 
1) 

Be epn sne lac tAgk Paglogey 
26 Z 


Cos wv -Sinw 0 


[AJ] = Sin w Cos wv 0 
0 0 1 
Cos 0 0 Sin 6 
[A.] = 0 i 0 
-Sin# 0 Cos 6 
sf 0 0 
[Ag] = 0 Cos ¢ -Sin 9 


0 Sin 9¢ Cos 6 


SUBROUTINE ROTAT calculates the transformation matrices for 
every NAMELIST PICT, except if ISCALE = 0. ISCALE = 0 


directs the scale of the plot to be the same as that of the 
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BEeev2ous Plot... Should.a_rotation. occur with ISCALE = 0, 
the plot width could exceed the paper width. 
be. scaling 

The safest scaling method is automatic scaling 
(ISCALE = 1). The user may specify a scale (ISCALE = 2) 
and the plot origin (XORGN,YORGN), but one must be careful 
not to run the plotting pen off the graph paper. ISCALE = 0 
is a very useful option. The plot will use the same scale 
as the previous plot. It is useful in an assembly graph 
where examination of a mesh in sections without losing 
perspective is desired. Example 3, Section III, illustrates 
the option ISCALE = 0. When ISCALE = 1 in a NAMELIST PICT 
defining a partial plot, a "blow-up" of that section is 
obtained. ISCALE cannot be zero in the first NAMELIST PICT. 

G. sPareral Plots 

To develop a partial plot, three methods of 
segregating elements exist: first, by the X, Y, and Z 
cutting planes; second, by node numbers, and, third, by 
element numbers. If a model has an area where the elements 
are relatively small, a "blow-up" may be desired. Choose a 
numbering scheme or coordinates to define the section to 
be segregated using one of the methods above. Example 3 
fsiguce £4, Section fii) uses xX, Y, and Z cutting planes to 
define the partial plots. Example 4 (figure 17, Section III) 


uses element numbers to section the plots. 
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C. NODAL POINT (GEOMETRY) INFORMATION READ-IN 

Nodal point data is read in by the GEOMn subroutines 
(GEOM1, GEOM2, and GEOM9, see figure 2). Since SAP IV 
and ADINA data decks are similar, SUBROUTINES GEOM1 and 
GEOM9 are also very similar. They are both constructed to 
read and generate data in exactly the same way as ADINA and 
SAP IV. All data not needed by PSAP1 is disregarded and 
the nodal point data-is stored in array 222 (figure 3). 
After studying storage array ZZZ and GEOM1 (or GEOM9), 
a user moderately familiar with FORTRAN programming could 
easily construct a user supplied subroutine (GEOM2) to 


read the nodal point data in any desired format. 


D. ELEMENT (CONNECTIVITY) INFORMATION READ-IN 

After reading and storing the nodal point data, the 
element data must be read. The GEOMn subroutine will read 
the element control card (NPAR, references 1 and 2). 
SUBROUTINE ELTYPE (figure 2) calls the proper element sub- 
routine to read the element data specified on the element 
control card. If several groups of elements are to be 
read, the process is repeated until all of the element 
groups have been read. Although the nodal point data is 
stored in array 22ZZ, the element connectivity is read and 


stored on device 10 (disk). 


E. DISPLACEMENT DATA READ-IN FOR DISPLACEMENT POSTPROCESSING 
Displacement data may be read in (figure 1) by SUBROUTINE 


DATA9 (KDATA=9) or SUBROUTINES DATA1 or DATA5 (user supplied, 
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KDATA=1 or 5). When read, displacement data is stored 

in the last three columns of array 22Z (figure 3). PSAP1 

can postprocess displacements for both ADINA and SAP IV. 

The difficulty comes in obtaining a punched deck of cards. 
Reference 6 gives a description of how to obtain a deck of 
cards for SAP IV in a format acceptable to SUBROUTINE DATA9. 
ADINA has no such provision. However, when preprocessing, 

the displacement data will be omitted (NUDISP=NVDISP=NWDISP=0) , 


and this step will be by-passed. 


F. PLOTTING LOGIC 

SUBROUTINE PLOTX (figure 4) is the main plotting routine. 
Since the nodal point data is stored in array 2ZZZ and the 
connectivity is stored on device 10, it is a simple matter 
to read the connectivity from device 10 (one element at a 
time), and connect the nodes as they are defined in references 
1 and 2. For example, NEND = number of nodes defining the 
connectivity of a single element, NUMEL = the element number 
and NODE(NEND) is the array containing the connectivity. 
Device 10 contains this information successively for each 


element. The statement 


READ (10) NEND,NUMEL, (NODE (I) , 1=1,NEND) 


will read the element connectivity to be plotted. MThe 
8-20 node brick (ADINA and SAP IV) and the 4-8 node plane 
elements (ADINA) may have 3 points defining each edge. If 


the midpoint node is defined, then isoparameteric 


24 


_ ; eS ) a es 


‘ oo 
betete ei stab Ineme ae , 


7. ove 
oie 7 site esl od? 
) - 
‘VI IAG bags AVIGA tidod, to3 sncmada nee aaeeeeee 
- 
; 


7 
Icsey .(£ enepid) 82S qadee 6 amie 


oe 
j j f i am t 
,sivtsoa Yo. > b berutve paiGists2o nf esiios ytivolt2ss 
Ss 


) feeb cs nisddo’ oy word So notdqtzoeet « wevig ay: 


a aurruo#rave [dstgs: dmetot » at VPSdAR- 10% 1 
on teses sic -fedw. .sevewoH .noletvo1wy Ac 
. ve | CWHi=fETQVAK Tet aur Beszsimo sti Iiliw #2s6 
= ce 
enitvos poig j, atem sat f ) ATO 
eas Oi: so YSRZIB cs ee 23 : 
-. TazJan Signse BSB ti 32 ,.9T : te 
& 2s inemeis ene) { sory mos. vos 
: SotetaIe: ai Bbent=eh. s2 atis ao 
eds patniieh eabor igduus = CMa ,elom 
atin a sir = ait , nae if@ ¢& 


° 
=. 
x 
. 
tt". 
‘. 
¢ 
' 
> 
re | 
‘3 
| 
~ 5 
‘4 
“ 
ba | 
5 
t 
Ww 
ry 
= 
i 
i 
5 
bs 


‘ 
mr 
Ait 
ny 
Yu 
Q 
{ 
4 
© 
ia) 
* 
is 
= 
- 
( 
La 
i) 
i 
‘ 
_ 


2nemecrete off Shem 


7 7 & 
(\dWav.f<2, (1).adom) .aauy cman (peat 


i 


wif .betsolg ae oo ‘ven +oec ote sé vs 


« 
yy 


w4 - 
anata shon 84+ wz) Bae (ys) ana] ine AUTGAY oboe 
7 j a es 
; 316 6.epbe Agse 1422386 mio £ eoetige ee. ‘“7GL) #30m 
= 7. er e& _ = , = J ? 
_— 7 : ee : — 


m a ee = 
; > hued bie tao aes: a3e: ot ebom smbeqbt 


START 


Read elements individually to be 
plotted from device 10. 
READ(10) NEND,NUMEL, (NODE(I) ,I=1,NEND) 


| Check if element is in bounds ne 
| to be plotted 


) YES 


NEND=2 NEND=4 NEND=8 NEND=12 OR 20 


Pot. 2 PLOE SAP TV Plot ADINA 


node 8 node 8 node plane 
elements elements brick and 20 node 
elements brick elements 
and SAP IV 
20 node brick 
elements 


NO Check if all elements have been 
read from device 10. 


YES 


RETURN 


Figure 4. Flow chart for PSAP1 plotting subroutine, PLOTX. 
Array NODE contains the connectivity of the 
element being plotted. NEND = the number of 
nodes defining the connectivity of a single 
element. NUMEL = the element number. 
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shape functions (reference 3) 


the three-node edges. 


to those used by ADINA and SAP IV, 


are used to interpolate along 


These shape functions are identical 


so the geometry represented 


graphically is identical to the problem solved in ADINA and 


SAP IV. 


following equations: 


Xo = Nl * Xa + N2 
Xo = Nl * You + N2 
we —-S F (S-= 1.0) 
ne =.—) (S_+ 1.0), * 
nS —S.* (5 fF 1.0) 
=i es ey Oe 


PSAP] uses the NPS 


subroutines. 


plotting package 
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III. PSAP1 SAMPLE PROBLEMS 


The following examples have been chosen to illustrate 
some of the most useful options of PSAP1l. Prior to 
attempting to use PSAP1, the user should have the problem 
defined and the cards punched in the format of references 
land 2. Appendix A of this thesis gives a complete des- 
cription of deck preparation for PSAPl here at NPS. This 
section should prove helpful in the understanding and 


interpretation of the options presented in Appendix A. 


A. ADINA EXAMPLES 

in heigrorced Concrete Beam (example 1, figure 5) 

This example was chosen because it illustrates 

the importance of the exploded plot when more than one 
element group is used. It is taken from the ADINA manual 
[reference 1]. Figure 6 is a listing of the data cards 
as they are prepared for ADINA. Figure 7 indicates how 
that deck would be modified for PSAP1. Note: load cards 
are removed, NAMELIST OPTION and NAMELIST PICT are added, 
and the title to be plotted on the graph is added in figure 
7. Otherwise, figures 6 and 7 are the same. Figure 8, 
part (a), illustrates an undistorted (KDISP = 0) PSAP1 plot 
with the nodes numbered (NOTAT = 1). Figure 8, part (b), 
shows the same mesh in an exploded form (KDISP = 2) with 
the elements numbered (NOTAT = 2). Note how the truss 


elements are visible in the exploded plot. 


2 


+4 


- anaseond, 


etszsanll t ‘oe me2ornS asad 

1679 AY ie ilies iweb a 
melddac sax seved- Sfoods tees “ana \iGhes 
esonsistes to) dentro ont nt bedsavd’ eb=d> oie 
-deh sJselamon-n aavip atzaas ei? 10 4 xiaenah’: a 
te.a3en 1 GARd 108 molt a upoSste toed tou sid 


A xibdeigé Ad. beenSaetq snokzao |eHe 
{© sot yf SlqmszS }, mage 2¢2s:506> beszoAnt ue 
= ete 2 Reese _leieaatieinellane ete aT ae a 


eter zev LLt sts esnebed: fSsons 2ew cicguisxe al 


v 
2 
Y 


si0om monw sole pebolgxue si? to sonsag 
Letineadg ARUASS AS mes uSAES ei st teas aF 
=buad seh) ert zo pafosig s ef os: Sireit «tL 
wa Badeod bret Sxip Ls ARICA ae2 Setegetg 
2bus5 beci. 73808 + f¢Aes no? ‘ekg foe od bigo 
Sebbs. exe TOTS rerianny’ Bete: § ; 


dolq lgaet (0 = cegauy 
(2) gaag .F eu it rr 
aigiw (f = sede 

vals siz wor asiht™ 


P72 P/2 


P/ 
wee 2 
2g A,, 20.62 inz P/2 


BEAM DIMENSIONS 


9 15 ays 
3 
r 
n 2 
4 
; 7 Is 19 25 3 
ke— 20'—eHe— 18" ee 2¢ oh Sele 


FINITE ELEMENT IDEALIZATION 


MATERIAL PROPERTIES: 


Ge = S40_psi 
co; = 458 psi Figure 5. Example 1, ADINA 


Z , truss and 8 node plane 
Py steel = 44000 psi elements, Reinforced Concrete 


3 = 6100 ksi Beam (Given on page 84, 
Sconcrete reference 2). 
y = O02 
Estee: 7 30000 ksi 
zy = 300 ksi 
steel 
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Pees INGt Wear Oi motte deck. 


‘ 
I 
i 
i 
i 
i 
i 
; i 
"MY PAP era rere ray 
4 a 
* 
* ete 
a Sonik ien @ 8 & & @ @ & 
i] 
uy ome fe - mein PE So LA Eee ol) = 
Leb 
} 
‘ 
{ I, 
- = 
th 
s! 
res R ** @* ¢.* © © @ + « © © Cw & © @ @ i 
, 
= Maha par ot on Lom = "de 
oy Pr i Ys es ab 0 a é A 65 
i 
- _ f 
at , - 
1 ive l 
- - -— 
“4 . 
~ 5; 
‘ = 7 Y 
” - in e 
‘ ~ ‘ ) - 
- = 
7) a 
how 
P 4 
x] b. . 
© 7 
4) ] f ” 
al = ao ~# 
i ml " ; 
Ls > > > 
I - ha bed > ee 
- + ~~ bo hem te tomy oy me featy Men. Fae fame Pa Get bey Pty “reel, pnts sey leis ae oe 
_ a = - 
4 \ os eel 
} oa > > 
i “ 
va! 7 \ 
i ~ ~~ 
*) “ , \. » He 
i = iiiadintinmnimion tooo hat a at TT 
tole = 
1 
ee — 7 7 
ad m% 
i i 4] = 
: « 
i ond 7 7 
7 ) = roe hee Eee Fe me hee eet fomte nee Rim, grnee heey heey Pees: Smeg oe te heey fi 
i 
# vw ~ I s ' 
~ nt 
; - 
= I “ 
aa] r s 


- | [26/9 POIaIOOO-AdIO IOS 


i fi i 


c tS) SS Oven She QOS OSB 


7 ae i 
* ~v - 7 ' 
; a) \ _ 
. - 
- — { Sy PE fo le hr ts ay en ro bs 
j - 
+d 


as 
tient lie iene 


ea Pie 


“sep Anan fOr ace jOh wae 


°Z JO zg obed ‘butyAst{T yoop yndut WNIdv ‘T eTdwexq °*9 sanbty 


S*O=" FI € ike 
G°ET O78 
OZ OFT Gog Gao CUS C°8 0° 
(0) “{S B°4 O°d (eS O°T 0°O ; ay 
Vivd QVO1 GQsIllddV Xeoak ok ok ok 
0) 62 0 O€ aE 92 LE EE 
°9 9 I! T 9 OT 
0) G 0 9 8 e € 6 
°9 9 T 1 9 9 
0) Bd @) 62 Te Gd 9¢ ce 
=9 9 T 1 9 G 
0) 4 0 S Ih T CA 8 
°9 9 T T 9 T 
0S°O S000 *0O cO0°- oy Lath 8S4°0O 
OGe@ SOLS 
VOITETCOU0- 1 
1 G 0 € 9 ra O T OT 
LAdNI LNS3W315 WANILNUD G2 3K ek ok 
TO) 1 "0 T Le B2 Ol 
"0 € 1 *c I 47 it T 
“OOE O° “OO000E 
aO. SIO ee LUO - SS 
T € O T OT il 


iNdNI LNSW31343 SSNYL Mek ek 


30 


_ > ae 7 ee?) 


= 


ac . OF 
<= 
—_ 
Pte 

> a mm -_ — 
£ Oo - 

iy ™ atta | 
= ry 


y «4 
ker) oF mk 


M BE PLOTTED ON GRAPH 
NFORCED CONCRETE BEAM (ADINA =X) 


1 
F 
NAMEL TST OPTION 


PSAPL ee 


KIBLER AE NONLINEAR SNALYSIS OF A RE= 
site ole ke 3K 


1K KOK AK ok 


XAMPLE 


uss 


TK OE KOK 3K 


NONLINEAR ANALYSIS OF A REFINFO 


Oe 


ONTROL CARDS 


0 
25 


SILOOLTY 


MAMMA MMMMMe MMMAMNMN 


fap) ceeeeee 
© &® © © 8&6 ORM aH Ree eS 


ttt cade 


eoreoeree ese et eevee eee @ 
TOODONDOCDODIODAIOOCWO AND 
AANMMND ANMMEMO NMMNMNO 


ee A ee he ee oe fo ee 


NODAL POINT DATA 


es ete Stat tat ace 


Fe Aan Roam Ree Me ce Noe Be Be Be ee oe oe hee Roe ee eh ke | 


FE EAE 2K OK 


el Ielselejer. leloalelolejelsiolea)\elelo) 


Sates edeseaet eee eset eet eS Ss See Ss 


SFE MNOIAANNOFOIONANONA OS 
Sam ANN BNNOS 


el 


t deck. 


inpu 


Werte eK NOte pate Ol 


CUNTROL CARD 


LOAD 
PSAP] anput deck listing, page 1 of 2. 


A 


HK A AS OK 


we 
Ae >< 


Example 1, 


20 


Figure 7. 


vy: 
i; i 
ir 
- 
: ii 
‘eate Ales Lae Leal 
phy 
I 
mq * 
2 - 
a 
ad 
7 i 
i) : 
=! 
vi a 
’ 
€ 
: 
E 
, 
- — 
— & — 
: . ~~ > 
é 
» aieils = 
aA “36° 00- IOS Qe eee "% 5 < 
Cd bene 
. Bi 
a 
i! : 
on > 
~7™ Se we ee es fe te -— es Ge dees fee bee © ree ~- 
’ ic 
i : - 2. 
4 . 2 
i u YY. = 
: 2 —_ ~~ 7 
~*~ = Amer mA m*> Oe ip aS hy <2) os a 7 = 
| ote HL Piso -_ Oy el ee. Pe 


) 


"2 FO ¢ ebea “burgsty, Syoep andut pavsd 


yoop yndut jo 3Aed JON xxxxy 


(HiO tid) “Gado VdxXa) SEoTd LST tSiwivin yo 


(AYNLINYLS GaWeyC330NN) LIId LSITISWYN ok ok XK ok 
Q3AQwW3SeY SI VLVG GVO1 Gal lddV ok ok ok ok 2k 

0 62 0 Of& ce 

paz 9 I 

0) G 0 9 8 

ve 9 if 

0 82 6) 62 Te 

ie) 9 T 

O 4 C S id 

wis) 9 T 
0S°O S000 °0 CUUl= 

T S 0 € 9 c 
JNdNI LNSWS3753 WNANILNOD dz OK OK OK 

so 0) I (0) nt 

°C € T OW 

HOOK 
“O)-€98'e 
I € 

SIN3SW513 SSNUL a He ke aK 2k 


SNOILIGNO)D IVILINI HX Re ak ok 


1 
7 
Ol 
a 
0) 


So Pea WANRGO 


coor 


‘{T oTduexy 


N 


N 


1s 
+ 
0 


"ZL oanbty 
IJN33 
*0*°6=Z7S101d 
*Z=1VLON 
*0=300» 
*c=dS 10» 
*2°O=S9VWG 
OM 3 
GN33 
*1T=300% 
"GL" 8=7S5 10 Id 
*T=LVLON 
®E=145 Ay 
*2=Z7YGHy 
PO ivals 
Lz Se 
9 OT 
e 6 
9 9 
92 Ze 
9 G 
c B 
9 I 
AS BSGbF°O 
60) SOO? 
LETEOOO 1 
T Ol erica 
82 OT 
T T 
4 *"OQO00E 
00” ae 
T Ole Tf 
T 


32 


: 
TUM! THIN224 MULMTTHOD OS 
0 | | 


a $ o 


: b2.6 F000 90: a 
= 4 


GQ2VUMae #1 ATID AA) AEP peaA eewes 
7a) WeTe OP VRLa Patt) TI fet gohan eirae 


(T0359 GAaGO IGE Fy) 7 4 te 1 JR an avazee 


toeh joqal to tyusq tak **ees 


Sto & spsqg ,palszetl d6ob smal PSABT .1 afqmexe .{ ssapld 


e 


IL 


7@ 30) p obed 


‘syudeib 4ynd4no [dwsd 


‘T oTdwexq 


°g oanbtg 


=LVLON) perequnu sopou ‘(Q=dSIGM) eAnqZoONnaAAS powsrOCJepup 


( 


e 


) 


c Jit | 


°Z JO ¢@ obed ‘penutju0D 6°g eanbtg 


"(Z=LWLON) perequnu sjuseuete / (z7=dSTay) 


30Td pepotTdxg 


(q) 


34 


a a = 


rc 
' aan —_ ’ 


——s Se ae ~ 1 
h 


2. Flat Plate With Hole (example 2, figure 9) 

This is a well known problem with which one can 
calculate the stress concentration on a hole in a plate 
under axial tension. Figure 10 is a listing of the PSAPI1 
deck set-up. The mesh is composed of ADINA variable 4-8 
node plane elements. Notice on figure ll, parts (a) and 
(b), how the interpolating shape functions round off the 
3-node edges. Part (a) has the nodes numbered (NOTAT = 1). 
Part (b) has the elements numbered (NOTAT = 2), and illus- 
trates the use of the symmetry option (KSYMXY=KSYMXZ=1). 
The symmetry option enables one to obtain a picture of the 
complete plate even though the model only consisted of a 


quarter plate with proper boundary conditions. 


B. SAP IV EXAMPLES 

1. SAP IV Truss Problem (example 3, figure 12) 

Figure 13 is a listing of the PSAP1 data deck. 

Figure 14 indicates how multiple partial plots can be used 
to obtain a better representation of the model. Part (a) 
of figure 14 is the complete model. Part (b) is the left 
mate 1 xxXMAx = 50"). and part (¢) is the right half 
tn — 50°) XeMAX = 1.05820"). Notice also that for 
Parts (b) and (c), ISCALE = 0, which means succeeding plots 
have the same scale as the first. Had ISCALE equaled l 
(blow-up), then the width of the half view would have been 
the same as the complete model. Figure 14 size is limited 


by the NAMELIST PICT variable PLOTSZ. 
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Figure 9. Example 2, flat plate with a hole in 
tension. ADINA 4-8 node plane element. 


Thickness =e neh 

Young*s modulus = 30.0 X 1p” psi 
Poisson's ratio = .333 

Total elements = 18 

Total nodes = 44 


Distributed Load = 


NS 


1000:..0' psi. 
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Example 2, PSAP1 input deck listing, page 2 of 3. 


Figure 10. 


**kEKK Not part Of input deck. 
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Figure ll. Example 2, PSAP1 output graphs, page 1 of 2. 
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(a) Actual mesh, nodes numbered. 
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(b) Symmetric representation, elements numbered. 
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NAMELIST PICT (LEFT HAND SIDE) 
NAMELIST PICT (RIGRT HAND SIDE) 


PSAP] input deck listing, page 2 of 2. 
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Example 3, 


Figure 13. 
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Figure 14. Example 3, PSAP1 output graphs. 


(c) Right hand 
side. 


(a) Complete structure. 


(ob) Lett hand 
side. 


_2. Cylindrical Bar With Spherical Hole (example 4, 
figure 15) 


This mesh could have several uses. Two might be 
to calculate stress concentrations if the bar is under 
axtal load orto calculate loading if the void is under 
pressure (i.e., dispersed nuclear fuel pellet). Figure 16 
is a listing of the PSAP1 data deck. Figure 17, part (a), 
is a representation of the complete model. Parts (b), 
7, 1a), and (e) of Eigure 17 are partial plots of the 
total structure using options of the undeformed structure 
(KDISP = 0) with node numbering (NOTAT = 1), and the 


exploded plot (KDISP = 2) with element numbering (NOTAT = 
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a two inch diameter spherical void on the center 
line under pressure. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


With a little study of PSAP1 and figure 2, one realizes 
how simple it would be to expand PSAP1 to include virtually 
an unlimited number of geometry data input formats. Simply 
supply the GEOMn subroutines, element subroutines and 
modify SUBROUTINE ELTYPE if necessary. Another interesting 
project would be to incorporate PSAP1 into ADINA and SAP IV 
so that a mesh plot could be obtained in the data check 
mode only. Reference 4 contains routines for plotting stress 
contours which could be incorporated. Hidden line logic 
is available. It would be difficult but could be incorporated. 
Hidden line plots are pretty, but they are not very practical 
when preprocessing because in preprocessing it is advantageous 


to see all of the nodes. 
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APPENDIX A 


PSAP1 USERS MANUAL 


I. .NAMELISTS 


It may be useful for the user to review NAMELISTs in any 
good FORTRAN manual. However, a short description of the 
NAMELIST input is given here. Only columns 2 through 80 
of each card will be read. The computer expects to find a 
special delimiter symbol in column 2 of the first card 
followed by the NAMELIST name. The delimiter symbol for 
the IBM 360-370 series is the ampersand & ; other processors 
use the dollar sign $. Following the first card comes the 
list of variables and their assigned values. Some important 
points about the variables are: 

A. Variables and their values need not appear in any 
special order. 

B. Predefined variables need not appear in the list 
Zeit 2S not necessary to change the value. 

C. A comma should follow each assigned value. 

D. It is recommended that each variable appear on a 
separate card to facilitate change. 

The last card following the NAMELIST variables contains, 
beginning in column 2, the delimiter sign followed by the 
word END. For example, &OPTION (read NAMELIST OPTION) 


contains 12 variables. All are initially assigned default 
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values. But, suppose three of those values (NNDEST=200, 
KGEOM=9, YSPACE=2.0) do not apply to a specific problem. 
The &OPTION would be read as follows: 

_&OPTION 

_NNDEST=300, 

_KGEOM=1, 


_YSPACE=5.0, 
_&END 


Blank space in column l. 
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LLiSjenENPUT PROCEDURES 


The following sequence of cards is necessary to use 


PSAP1. 


NOTES IBM CARDS 
COLUMNS 


123456789-------------- 


(1) jf shoeandard green job card] 
// EXEC FORTCLGP, REGION.GO=150K 


F/EORESSYSIN DD * 


(2) E MAIN PROGRAM 
DIMENSION ZZZ(NZ) ,DISPD(5,3,NON) 
CALL PSAP1(ZZZ,NZ,DISPD,NON) 
STOP 


END 


(3) ae 
//LINK.USDD DD UNIT=3330,VOL=SER=DISKO02, 
// DISP=SHR,DSN=S1153.PSAP1 
/ / GaNigesvsin DD * 
INCLUDE USDD (LOADM) 


ENTRY MAIN 


62 


ie ee ty a 
eu Od vteseSoen ex ese Oo soaeups 
a = 


(buse dot sess bes bist 
= £ — co 

5. SEITAO Dax * 
J: ee 7 
yy 


: aa wteve.' LOT 


uve, 2 ORRIN RE worausmra 2 


= (Wi OW, age FG, eu, sen) giaee aad : 
. | i - j 


*< 4 gore” 4. 


(4) 


2) 


(6) 


(9) 


(10) 


(11) 


(12) 


/* 
//GO.FT1LOFO01 DD UNIT=SYSDA, 
Jf MSPACE= (CYL, (3,1)), 


// DCB= (RECFM=VS , BULKSIZE=3520) 


7/GO.SYSIN DD * 


[PSAP1 title card as it appears on plot] 


&OPTION 


[&OPTION variables to be initially set or changed] 


& END 


[ADINA or SAP IV geometry data. Title 
card through element data — remove load cards] 


[Case identification card — Omit if IDCASE = 0] 


[Displacement data cards — Omit for preprocessing 
only] 


&PiCcT 
[&PICT variables to be initially set or changed] 


& END 


[Additional data as defined by last value of 
KODE in NAMELIST PICT — Omit if last 
value of KODE = 0] 


/* 
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Notes: 


(1) Standard basic deck set up as described in chapter 
3 of reference 8. 


// (Standard green job card] 
// EXEC FORTCLGP,REGION.GO=150K 
A fEORT.SYSIN DD * 


(2) Main program. 
DIMENSION 22ZZ(NZ) ,DISPD(5,3,NON) 
CALL PSAP1(Z2ZZ,NZ,DISPD,NON) 
STOP 
END 
The main program has two functions: to allocate 
fast storage space, and to call PSAP1. NON must 
be greater than the number of nodes. NZ must be 
+3 
greater than 4*NON(7*NON if displacement data cards 


anew to be input, it.e., NUDISP, NVDISP, or NWDISP = 1). 


(3) Using load module library. 
/* 
//LINK.USDD DD UNIT=3330,VOL=SER=DISKO02, 
// DISP=SHR, DSN=S1153.PSAP1 
//GENK=SYSIN DD * 
INCLUDE USDD (LOADM) 
ENTRY MAIN 
PSAP1 should be precompiled and stored in the machine 
in a load module. Load modules are a type of user 
library described in reference 7 and chapter 3 of 


reference 8. Since PSAP1 requires over one minute 


to compile, precompiling results in a large time 
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(4) 


(5) 


Savings. These cards are subject to change, and the 
most current version is determined by the user main- 
taining this library. If PSAP1 is’ not on a load 
module, these control cards may be replaced by 
subroutine PSAP1 (which consists of a box and a 


half of cards). 


Allocation of storage. 
/* 
//GO.FTLOFOO1 DD UNIT<=SYSDA, 
V7 SPACE=(CYL, (B,1)%) , 
// DCB=(RECFM=VS , BULKSIZE=3520) 
PSAP1 uses a Slow storage device to store the element 


connectivity. These cards allocate 3 cylinders as 


described in chapter 3 of reference 8. 
Deck set up card. 
7 /GO.SYSIN DD ~* 
This is a standard deck set up card described in 
chapter 3 of reference 8. 


PSAPL title card. 


PSAPL title card is the title that will appear on the 
graph. Make sure a user identification is on this 
Gard. it consists of 80 alphanumeric characters. 

The first 40 characters will form the first title 


line. The last 40 will form the second line. 
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(7) NAMELIST OPTION. 


(Note: start in second column) 


&OPTION 
[SOPTION variables to be initially set or changed] 
& END 
VARIABLE-DEFAULT DESCRIPTION 
NNDEST-200 Must be equal to the number of nodes. 
NUDISP-0 0 — X direction displacements not input. 
1 — X direction displacements input. 
NVDISP-0 0 — Y direction displacements not input. 
1 — Y direction displacements input. 
NWDISP-0 0 — Z direction displacements not input. 


1 — Z direction displacements input. 
(Note: unless displacement data is to be input, 
allow NUDISP, NVDISP, and NWDISP to default.) 
KGEOM-9 Specifies the geometry input format. 
1l- Subroutine GEOM1 reads in ADINA data 
deck geometry and connectivity. 
2- Subroutine GEOM2 may be supplied by 
user along with subroutines NSPLAN, NS3DEE and NSTRUS 
to read nodal data and connectivity in any format. 
9- Subroutine GEOM9 reads in SAP IV 
data deck geometry and connectivity. 
KDATA-9 Specifies the subroutine and corresponding 
method of input for displacement data. 
1- Subroutine DATAl1, a user supplied 


subroutine. 
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2- Subroutine DATA2, a user supplied 
subroutine. 

9- Subroutine DATA9, reads a punched 
output displacement deck from execution of SAP IV as 
presented in reference 6. (Note: unless displacement 


data is to be plotted, allow to default.) 


NVALUS -0 Not incorporated, allow to default. 
IRESEQ-1 Not incorporated, allow to default. 
KPLOT-1 Not incorporated, allow to default. 
YSPACE~-2.0 Space between plots in the Y direction 


in inches when successive plots are plotted (i.e., 
KODE # 0). The graph title is plotted both on the 

top and at the bottom of each set of graphs controlled 
by a given NAMELIST OPTION. The space between the 
elebe and the plot is YSPACE/2.0. 

PSIZE-9.0 Paper size in the X direction, in 
inches. Used in scaling of the plots to insure this 
dimension is not exceeded. However, when manual 
scaling (ISCALE = 2; see NAMELIST PICT) this protection 
is not available, and it is possible to exceed the 
paper width. 

IDCASE-0 O- No identification card precedes 
the deck of displacement values. 

l- Identification card precedes the 


deck of displacement values. 
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SAP IV or ADINA data cards. 


Here insert the geometry (node coordinates) and 
element connectivity. This includes the title card 
through the element data cards; the load cards are 
removed. For a SAP IV data deck, NAMELIST OPTION 
variable KGEOM = 9. For an ADINA data deck, 


KGEOM= 1. 


From ADINA deck remove: 
l. Applied loads data. 


2. Frequency and mode shape calculations data. 


From SAP IV deck remove: 
1. Concentrated load mass data. 
2. Element load multipliers. 


3. Dynamic analysis cards. 


Otherwise, these cards are exactly the same as the 
deck prepared for SAP IV or ADINA. PSAP1 is not 

limited to SAP IV and ADINA. The user may specify 
any unique format (i.e., KGEOM = 2). In this case, 


SUBROUTINE GEOM2 must be supplied by the user. 


Case ID card and displacement data cards. 


PSAP1 is intended to be used essentially for prepro- 


cessing. If it is desired to use the postprocessing 


option, the user is referred toreference 6. Otherwise, 


omit the case ID card and displacement data. 
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Pre) and (11) NAMELIST PICT. 


(Note: start in second column) 


&PICT 
[&PICT variables to be initially set or changed] 
&END 
VARIABLE-DEFAULT DESCRIPTION 
KHORZ-1 Integer designating the horizontal axes 
of the viewing plane. 1 = Xo: 2 = Yo: 3 = ae 
(See figure 18.) 
KVERT-2 Integer designating the vertical axes 
of the viewing plane. 1 = X5- 2 = be 3 = oe 
(See figure 18.) 
PHI-0.0 Angular rotation Gf the model@about its 


X axis in degrees (performed 3rd, see figure 18). 
THETA-0.0 Angular rotation of the model about its 
Y axis in degrees (performed 2nd, see figure 18). 
PSI-0.0 Angular rotation of the model about its 
Z axis in degrees (performed lst, see figure 18). 
NEWFR-1 l- Frame change before plotting. 
2- No frame change before plotting. 
(Normally allow to default. A frame change resets 
the Y origin past the previous plot by YSPACE given 
in NAMELIST OPTION and resets the X origin at 0.0) 
ISCALE-1 O- No scale change. Use the same scale 
as the previous plot. Useful in an assembly graph 


where it is desired to examine a mesh in sections 
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without losing perspective. See example 3 in 
Section III of this thesis. ISCALE cannot be zero 
in the first NAMELIST PICT. 


i= Automatic scaling of plot and 


computation of proper origin location. 


2- User specified origin and scaling. 


| PLOTSZ-10..0 Maximum dimension desired on com- 
pleted plot prior to rotation. After rotation it 
is possible for the maximum dimension in the pro- 
jected plane to exceed PLOTSZ. The maximum size 
of the projection is unlimited on the vertical axis 
but is limited by the paper width (PSIZE in NAMELIST 
OPTION) on Che horizontal axis. If rotation of the 
model causes the projection to go off the paper, 
it is rescaled prior to plotting. PLOTSZ is used 


im scaling only 1£f ISCALE = 1). 


XORGN-0.0 ~ location of the plot origin. 


(Used only if ISCALE = 2.) 


YORGN-0.0 WVelocaeion of the plot origin: 


(Used only if ISCALE = 2.) 


PSCALE-1.0 Model size reduction factor. 
PSCALE is equal to the actual model size divided by 
the desired plot size. It is used only if ISCALE = 2. 
(Note that when ISCALE = 2 is used, no rescaling 
occurs if rotation causes the projection to exceed 


the paper width.) 
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NOTAT-0 O- No numbering on plots. 
1- Numbering of grid points. 


2- Numbering of elements. 


XLHT-0.14 Height in inches of the integers 
specified by NOTAT. It should be a multiple of 
Wg at XL US, not a. muLeiple of ..07, 26 will be 
rounded to the nearest multiple. XLHT has a 


Maximum of .49 and a minimum of .07. 


KDISP-0 O- Plot of undeformed structure. 

im PLOt Of Geftormed SeEructure. 

2- Exploded plot. 

3- Displacement represented by 
VECTORS. 

KDISP = 1 or 3 represents a form 
of postprocessing and displacement data must be 
input in (9). Lf postprocessing) is) desined,’ reter 


to reference 6. 


IDMAG-2 1- Direct magnification of displace- 
ment data by DMAGS. 
2- Scaling of displacement data 


to a maximum value of DMAGS. 


DMAGS-1.0 Magnification of displacements 
(feeepESe = Lor 3). Reduction factor of elements 
Gis KDISP = 2). 

KSYMXY-0 1- Symmetry about X-Y plane. 


i2 


»eoa:oe & 

°* 
= ” ” 
- + ~~ 


siacizium 6 Toa & 
. ay 7 


s eed The -sigeslar tyetaen edt oF 


o. 20 aumicio & fine @6.°30 git 


> 


—— £. » 
15 5 tele 36 7olg =! ; 
- A : = = at ~Ol Gx ~*. 

A : 

: / 2 3 s96ice | @=€ he 

al 
a A 
7 ys @7070ey 
i e Al 7 . 7 = 

j aol 5 ~ 3 _— - 4 — - — err \ - +2 
; of = avast os £ ao { = @eray : 


a © m 6ta Jmomaneiaet® cas eainaesoryseeg to 


1eie7 . bexias : oleescsosgtaog 32 (2) a2 9 nk 
= - 7 
: av aousze3 ox & 


iSinpen 2soseii0 


ssloels toa pratisead 


» PQAMKT te. hg rat cuba ao 


é + 
ee 0958, Fale Ie wo fents ny : i: — 
eae : va i : 


_ 


* 2hemes T- “asd nets a hs i = atone 
piss os 


- ¢ : . " ue #ae2aN 


i « 


! - : ie 


2 ‘ef ag ¥-Z Soha : y 
oe: 


an 


KSYMXZ-0 l- Symmetry about X-Z plane. 


KSYMYZ-0 l- Symmetry about Y-Z plane. 
A plate quadrant with KSYMXZ and 
KSYMXY equal to 1 would yield a complete plate. 


See example 2 in Section III in this thesis. 


Note: To develop a partial plot, three methods of 
segregating elements exist: First, by X, Y and Z 
cutting planes; second, by node numbers; and third, 
by element numbers. The next ten variables are used 


to separate elements into partial plots. 


XXMAX , YYMAX, ZZMAX-1.0E20 Local cutting planes 

Pein rE AZ BEATIN Wb, OO) ee et FO the principal 
planes. 

NDMAX-9999999 Maximum gridpoint identification 


number to be included in the plot. 


NDMIN-O Minimum gridpoint identification 


number to be included in the plot. 


NELMAX-9999999 Maximum element number to be 


ine Lluded nm the ‘plot . 


NELMIN-O Minimum element number to be 


included in the plot. 
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KODE-0 


Specifies the control option after 
the plot is completed. 

O= Last*plot, exit from program. 

l- Read another NAMELIST PICT. 

2- Read a new set of displacement 
data. (Postprocessing only.) For KODE = 2, dis- 
placement data must be followed by another NAMELIST 
Pre: 

3- Read a complete new set of input 
data Starting with a title card. 

For KODE = 1, 2 or 3, addreional 
sections of the deck must be repeated. The deck 
must end with a NAMELIST PICT having a value of 


KODE = 0. 


Note: A most important point to remember when 
generating a sequence of plots is that once a 
parameter has been defined, it retains that value 
until it is reassigned. For example, if PLOTSZ is 
assigned a value of 8.0 for the first of a series 
of plots, and it is not redefined in any subsequent 
NAMELIST PICT; the value of PLOTSZ will be retained 
as originally specified. However, when KODE = 3 
and a new title, NAMELIST OPTION and problem data 
are read, all variables in NAMELIST OPTION and 


NAMELIST PICT are assigned their default values. The 
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problem starts over in this case. Refer to the 
flow +chart in figure il sand study the path for 


different values of KODE. 


Delimiter card. 


/* 


~ 


Delimiter card is defined in chapter 3 of reference 
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III. SPECIAL FEATURES OF PSAP1 


A. POSTPROCESSING 

Reference 6 contains information and examples on the 
use of the postprocessor. SAP IV has the capability to 
punch displacement data cards in an acceptable format for 
PSAP1 (specifically subroutine DATA9). As of this writing, 
ADINA does not have this capability, and data cards would 


have to be punched manually by the user. 


B. PARTIAL DATA 

Reference 6, on page 30, establishes a procedure by 
which it is possible to input only a portion of the finite 
element model for a data check. This feature is valuable 
in a case where several different people are preparing 
different parts of a large data base for a problem and 
desire to check individually their inputs graphically for 
accuracy. All nodal coordinates for the entire model may 
be input, or only those that specifically define the portion 
of the finite element model to be plotted. In either case, 
all nodal coordinates that define the elements to be plotted 
must be specified. To use this feature the "element control 
cards" (described in references 1 and 2) must be modified. 
For example, if only the connectivity for elements 15 
through 50 of a problem are available, the changes below 


would be made. 
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1. All SAP IV Elements 
All SAP IV element control card changes would 
have similar changes because columns 6-10 contain the 
total number of group elements. Columns 66-70 are not 
used. Thus to plot only elements 15 through 50, make the 
following changes: 
a. Enter the upper bound (i.e., 50) in columns 
6-10. 
b. Enter the Lower “bound” (i.e.7 15) im ecolumns 
66-70. 
2. ADINA Truss, 2D and 3D Elements 
ADINA element control cards for the truss, 2D 
continuum and the 3D continuum elements would be changed 
as follows: 


as hiecer thes upper bound’ (2 .e.7 50) “in columns 


b. Enter the lower bound (i.e., 15) in columns 
30:6 
3. ADINA Beam Element 
ADINA element control cards for the beam element 
would be changed as follows: 


a. Enter the upper bound (i.e.,, 50) in columns 


me enter the lower bound>(i.e., 15) in columns 


65-69. 
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APPENDIX B 
SUBROUTINE PSAP1 LISTING 
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